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ABSTRACT 

Based on an analysis of what it may mean for one tensor to depend in the proper way 
on another, I prove that, under certain natural conditions, there can be no tensor whose 
interpretation could be that it represents gravitational stress-energy in general relativity. 
It follows that gravitational energy, such as it is in general relativity, is necessarily non- 
local. Along the way, I prove a result of some interest in own right about the structure 
of the associated jet bundles of the bundle of Lorentz metrics over spacetime. 
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1 Gravitational Energy in General Relativity 

There seems to be in general relativity no satisfactory, localized representation of a quantity whose 
natural interpretation would be "gravitational (stress-)energy" . The only physically unquestionable 
expressions of energetic quantities associated solely with the "gravitational field" we know of in 
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general relativity are quantities derived by integration over non-trivial volumes in spacetimes satis- 
fying any of a number of special conditions; 1 These quantities, moreover, tend to be non-tensorial 
in character. In other words, these are strictly non-local quantities, in the precise sense that they 
are not represented by fields defined at individual spacetime points (such as tensors or scalars). 2 

This puzzle about the character and status of gravitational energy emerged simultaneously with 
the discovery of the theory itself. 3 The problems raised by the seeming non-localizability of gravita- 
tional energy had a profound, immediate effect on subsequent research. For example, it was directly 
responsible for Hilbert's request to Noether that she investigate conservation laws in a quite general 
setting, the work that led to her famous results relating symmetries and conservation laws. 

Almost all discussions of gravitational energy in general relativity however, dating back to those 
earliest debates, have been plagued by vagueness and lack of precision. The main result of this paper 
addresses the issue head-on in a precise and rigorous way. Based on an analysis of what it may mean 
for one tensor to depend in the proper way on another, I prove that, under certain natural conditions, 
there can be no tensor whose interpretation could be that it represents gravitational stress-energy in 
general relativity. It follows that gravitational energy, such as it is in general relativity, is necessarily 
non-local. Along the way I prove a result of some interest in own right about the structure of the 
associated first two jet bundles of the bundle of Lorentz metrics over spacetime. 

2 The Principle of Equivalence: A Bad Argument 

The most popular heuristic argument used to attempt to show that gravitational energy such as 
it is, either does not exist at all or does exist but cannot be localized invokes the "Principle of 
Equivalence". Choquet-Bruhat (1983, p. 399), for example, puts the argument like this, 

This 'non local' character of gravitational energy is in fact obvious from a formulation of 
the equivalence principle which says that the gravitational field appears as non existent 
to one observer in free fall. It is, mathematically, a consequence of the fact that the 
pseudo-riemannian connexion which represents the gravitational field can always be made 
to vanish along a given curve by a change of coordinates. 

Goldberg (1980, pp. 469-70) makes almost exactly the same argument, though he draws the conclu- 
sion in a slightly more explicit fashion: 5 

1 Weyl (1921, pp. 271-272) was perhaps the first to grasp this point with real clarity. See also Dirac (1962). 

2 Schr6dinger (1950, pp. 104-105) gives a particularly clear, concise statement of the relation between the fact 
that the known energetic, gravitational quantities are non-tensorial and the fact that integration over them can be 
expected to yield integral conservation laws only under restricted conditions. 

3 The first pseudo-tensorial entity proposed to represent gravitational stress-energy dates back to Einstein (1915), 
the paper in which he first proposed the final form of the theory. 

4 Sce, e.g., Brading (2005). 

5 Goldberg's formulation of the argument exhibits a feature common in the many instances of it 1 have found in the 
literature, the conclusion that a local gravitational energy scalar density does not exist and not that a gravitational 
stress-energy tensor does not exist. Perhaps one could imagine having a well-defined scalar energy density of a field 
in the absence of a well-defined stress-energy tensor for that field, though 1 cannot myself see any way to represent 
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[I]n Minkowski space any meaningful energy density should be zero. But a general space- 
time can be made to appear Minkowskian along an arbitrary geodesic. As a result, any 
nontcnsorial 'energy density' can be made to be zero along an arbitrary geodesic and, 
therefore, has no invariant meaning. 

Trautman (1976, pp. 135-6) has also made essentially the same argument. In fact, the making of 
this argument seems to be something of a shared mannerism among physicists who discuss energy 
in general relativity; it is difficult to find an article on the topic in which it is not at least alluded 
to. 6 

There are two fundamental flaws with arguments of this sort. First, the Principle of Equivalence 
does not deserve its definite article — there are almost as many formulations of it as there are people 
who invoke it, and many of those formulations are so vague as to leave it unclear whether or not 
the proposition holds in general relativity. 7 The second is that it assumes that, if there is such 
a thing as localized gravitational energy or stress-energy, it can depend only on "first derivatives 
of the metric". But that seems wrong on the face of it. If there is such a thing as a localized 
gravitational energetic quantity, then surely it depends on the curvature of spacetime and not on 
the affine connection (or, more precisely, it depends on the affine connection at least in so far as 
it depends on the curvature), for any energy one can envision transferring from the "gravitational 
field" to another type of system in a different form (e.g., as heat or a spray of fundamental particles) 
in general relativity must at bottom be based on geodesic deviation, 8 and so must be determined 
by the value of the Riemann tensor at a point, not by the value of the affine connection at a single 
point or even along a curve. There is no solution to the Einstein field-equation that corresponds in 
any natural way to the intuitive Newtonian idea of a constant gravitational field, i.e., one without 
geodesic deviation; that, however, would be the only sort of field that one could envision even being 
tempted to ascribe gravitational energy to in the absence of geodesic deviation, and that attribution 
is problematic even in Newtonian theory. 

An obvious criticism of my response to the standard line, related to another popular argument 
given for the non-existence or non-locality of gravitational energetic quantities, is that it would 
make gravitational stress-energy depend on second-order partial derivatives of the field potential 
(the metric, so comprehended by analogy with the potential in Newtonian theory), whereas all other 
known forms of stress-energy depend only on terms quadratic in the first partial derivatives of the 
field potential. To be more precise, the argument almost always runs like this: 

One can make precise the sense in which Newtonian gravitational theory is the 'weak- 
field' limit of general relativity. 9 In this limit, it is clear that the metric field plays roughly 

such an idea in general relativity. (Note that if one could, this would appear to be a violation of the thermodynamic 
principle that all energy is equivalent in character, in the sense that any one form can always in principle be tranformed 
into any other form). 

6 Bondi (1962) and Geroch (1973) are notable exceptions. I take their discussions as models of how one should 
discuss energetic phenomena in the presence of gravitational fields. 

7 See, for instance, Dicke (1962) and Norton (1985) for discussion of these points. 

s Pcnrosc (1966) and Ashtckar and Penrose (1990) implicitly rely on the same idea to very fruitful effect. 
9 See, e.g., Malament (1986). 
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the role in general relativity that the scalar potential 4> does in Newtonian theory. In 
Newtonian theory, bracketing certain technical questions about boundary conditions, 
there is a more or less well-defined energy density of the gravitational field, proportional 
to (V0) 2 . One might expect, therefore, based on some sort of continuity argument, or 
just on the strength of the analogy itself, that any local representation of gravitational 
energy in general relativity ought to be a "quadratic function of the first partials of 
the metric". 10 No invariant quantity at a point can be constructed using only the first 
partials of the metric, however, so there can be no scalar or tensorial representation of 
gravitational energy in general relativity. 

(No writer I know makes the argument exactly in this form; it is just the clearest, most concise 
version I can come up with myself.) As Pauli (1921, p. 178) forcefully argued, however, there can be 
no physical argument against the possibility that gravitational energy depends on second derivatives 
of the metric. Just because the energy of all other known fields have the same form in no way implies 
that a localized gravitational energy in general relativity, if there is such a thing, ought to have that 
form as well. Gravity is too different a field from others for such a bare assertion to carry any 
weight. As I show in §6, moreover, a proper understanding of tensorial concomitants reveals that 
an expression linear in second partial derivatives is in the event equivalent in the relevant sense to 
one quadratic in first order partials. This illustrates how misleading the analogy with Newtonian 
gravity can be. 

3 Geometric Fiber Bundles and Concomitants 

The introduction of a coordinate system to geometry is an act of violence. 

Hermann Weyl 

Philosophy of Mathematics and Natural Science 

I have argued that, if there is an object that deserves to be thought of as the representation of 
gravitational stress-energy of the gravitational field in general relativity, then it ought to depend 
on the Riemann curvature tensor in an appropriate sense. Since there is no obvious mathematical 
sense in which a general mathematical structure can "depend" on a tensor, the first task is to say 
what exactly this could mean. I will call a mathematical structure on a manifold that depends in 
the appropriate fashion on another structure on the manifold, or set of others, a concomitant of it 
(or them). 

As near as I can make out, this terminology is due to Schouten. 11 The definition Schouten 
proposed — the only one I know of in the literature — is expressed in terms of coordinates: depending 

10 In this light, it is interesting to note that gravitational energy pseudo-tensors tend to be quadratic in the first-order 
partials of the metric. 

11 See Schouten (1954, p. 15), though of course he used the German Komitant. The idea of proving the uniqueness 
of a tensor that "depends" on another tensor, and satisfies a few collateral conditions, dates back at least to Weyl 
(1921, pp. 315-18) and Cartan (1922). In fact, Weyl proved that, in any spacetime, the only two-index symmetric 
covariant tensors one can construct at a point, using only algebraic combinations of the components of the metric and 
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on what sort of concomitant one was dealing with, the components of the object had to satisfy 
various conditions of covariance under certain classes of coordinate transformations. This makes it 
not only unwieldy in practice and inelegant, but, more important, it makes it difficult to discern 
what of real significance is encoded in the relation of being a concomitant in particular cases. 
Schouten's invariance conditions translate into a set of partial differential equations in a particular 
coordinate system, which even in fairly trivial cases turn out to be forbiddingly complicated. 12 ft 
is almost impossible to determine anything of the general properties of the set of a particular kind 
of concomitant of a particular object by looking at these equations. I suspect that it is because 
these conditions are so complex, difficult and opaque that use is very rarely made of concomitants in 
arguments about spacetime structure in general relativity. This is a shame, for the idea is, I think, 
potentially rich, and so calls out for an invariant formulation. 

The reason I am inquiring into the possibility of a concomitant in the first place, when the 
question is the possible existence of a representation of gravitational stress-energy tensor, is a simple 
one. What is wanted is an expression for gravitational energy that does not depend for its formulation 
on the particulars of the spacetime, just as the expression for the kinetic energy of a particle in 
classical physics does not depend on the particular interactions one imagines the particle to be 
experiencing with its environment, and just as the stress-energy tensor for a Maxwell field can be 
calculated in any spacetime in which there is a Maxwell field, irrespective of the particulars of the 
spacetime, in contradistinction to the definitions of all known expressions for gravitational energy 
in general relativity now do (e.g., the ADM mass, which can be defined only in asymptotically flat 
spacetimcs). If there is a well- formed expression for gravitational stress-energy, then one should 
be able in principle to calculate it whenever there are gravitational phenomena, which is to say, 
in any spacetime whatsoever — it should be a function of the curvature in that spacetime, in some 
appropriately generalized sense of 'function'. This idea is what a concomitant is supposed to capture. 

its first two partial derivatives in a coordinate system at that point, that are at most linear in the second derivatives of 
the metric, are linear combinations of the Ricci curvature tensor, the scalar curvature times the metric and the metric 
itself. In particular, the only such divergence-free tensors one can construct at a point are linear combinations of 
the Einstein tensor and the metric with constant coefficients. Using Schouten's definition of a concomitant, Lovelock 
(1972) proved the following theorem: 

Let (JVC, g ab ) be a spacetime. In a coordinate neighborhood of a point p 6 JVt, let a( 3 be the components 
of a tensor concomitant of {g\^; gx^y, gx^^p} such that 

v n e nb = o. 

Then 

©ab = rG ab + qg ab , 

where G ab is the Einstein tensor and q and r are constants. 

This is a much stronger result in several ways than Weyl and Cartan had been able to attain: one has a more 
generalized notion of concomitant than algebraic combination of coordinate components; one does not demand that 
© a 5 be symmetric; and most strikingly, one does not demand that Q ab be at most linear in the second-order partial 
derivatives of the metric components. 

12 For a good example of just how hairy these conditions can be, sec du Plessis (1969, p. 350) for a complete set 
written out explicitly in the case of two covariant-indcx tensorial second-order differential concomitants of a spacetime 
metric. 
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The machinery of fiber bundles is used to characterize the idea of a concomitant. I give a (brief) 
explicit formulation of the machinery, because the one I rely on is non-standard. (We assume from 
hereon that all relevant structures, mappings, etc., are smooth. Nothing is lost by the assumption 
and it will simplify exposition. All constructions and proofs can easily be generalized to the case of 
topological spaces and continuous structures.) 

Definition 3.1 A fiber bundle 03 is an ordered triplet, 03 = (23, M, ir), such that: 
FBI. 23 is a differential manifold 
FB2. M is a differential manifold 
FB3. 7r : 23 — > M is smooth and onto 

FB4. For every g,p£ M, it~ l {q) is diffeomorphic to 7r _1 (p) (as submanifolds ofT>) 

FB5. 23 has a locally trivial product structure, in the sense that for each q E M there is a 
neighborhood U 3 q and a diffeomorphism £ : 7r _1 [[/] — > U x ir~ l {q) such that the action of n 
commutes with the action of £ followed by projection on the first factor. 

3 is the bundle space, M the base space, tt the projection and 7r _1 (g) the fiber over q. By a convenient, 
conventional abuse of terminology, I will sometimes call 23 itself 'the fiber bundle' (or 'the bundle' 
for short). A cross-section k is a continuous map from M into 23 such that ir(n(q)) = q, for all q in 
M. 

In particular, this definition of a fiber bundle is non-standard in so far as no group action on the 
fibers is fixed from the start; this implies in particular that no correlation between diffeomorphisms 
of the base space and diffeomorphisms of the bundle space is fixed. 13 On the view I advocate, 
the geometric character of the objects represented by the bundle arises arises not from the group 
action directly, but only after one fixes an explicitly defined correlation between diffeomorphisms 
on the base space with those on the bundle space. For example, depending on how decides that a 
diffeomorphism on the base space ought to induce a diffeomorphism on the bundle over it whose fibers 
consist of 1-dimcnsional vector spaces, one will ascribe to the objects of the bundle the character 
either of ordinary scalars or of n- forms (where n is the dimension of the base space). The idea is 
that the diffeomorphisms induced on the bundle space then implicitly define the group action on the 
fibers appropriate for the required sort of object. 14 

I call an appropriate mapping of diffeomorphisms on the base space to those on the bundle space 
an induction. (I give a precise definition in a moment.) In this scheme, therefore, the induction 
comes first conceptually, and the relation between diffeomorphisms on the base space and those 
they induce on the bundle serves to fix the fibers as spaces of geometric objects, viz., those whose 



13 In this I follow Geroch (1996). See, e.g., Steenrod (1951) for the traditional definition and the way that a fixed 
group action on the fibers induces a correlation between diffeomorphisms on the bundle space and those on the base 
space. 

14 1 will not work out the details of how this comes about here, as they are not needed for the arguments of the 



paper; see Curiel (2009). 
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transformative properties are tied directly and intimately to those of the ambient base space. This 
way of thinking of fiber bundles is perhaps not well suited to the traditional task of classifying 
bundles, but it turns out to be just the thing on which to base a perspicuous and useful definition 
of concomitant. Although a diffeomorphism on a base space will naturally induce a unique one on 
certain types of fiber bundles over it, such as tensor bundles, in general it will not. There is not 
known, for instance, any natural way to single out a map of diffcomorphisms of the base space into 
those of a bundle over it whose fibers consist only of spinorial objects. 15 

I turn now to making this intuitive discussion more precise. A diffeomorphism $ of a bundle 
space 23 is consistent with <j>, a diffeomorphism of the base space JVC, if, for all u e 23, 

7r(^(u)) = <j>(ir(u)) 

For a general bundle, there will be scads of diffcomorphisms consistent with a given diffeomorphism 
on the base space. A way is needed to fix a unique consistent with a so that a few obvious 
conditions are met. For example, the identity diffeomorphism on JVC ought to pick out the identity 
diffeomorphism on 23. More generally, if <f> is a diffeomorphism on JVC that is the identity inside an 
open set OcM and differs from the identity outside O, it ought to be the case that the mapping 
picks out a <jfi that is the identity on 7r _1 [0]. If this holds, we say that that <jr is strongly consistent 
with <j>. 

Let Dm and be, respectively, the groups of diffcomorphisms on JVC and 23 to themselves, 
respectively. Define the set 

®s = W *= ®2 • ^ *= ®M sucn that $ is strongly consistent with (/>} 

It is simple to show that Dig forms a subgroup of 25 £• This suggests 

Definition 3.2 An induction is an injective homomorphism i : X>m 

4> will be said to induce $ (under i) if i(<j>) = (jfi . 16 

Definition 3.3 A geometric fiber bundle is an ordered quadruplet (23, JVC, ir, t) satisfying 

GFB1. (23, M, tt) satisfies FB1-FB5 

GFB2. i is an induction 

Geometric fiber bundles are the appropriate spaces to serve as the domains and ranges of concomitant 
mappings. 

Most of the fiber bundles one works with in physics are geometric fiber bundles. A tensor bundle 
23, for example, is a fiber bundle over a manifold JVC each of whose fibers is diffeomorphic to the 
vector space of tensors of a particular index structure over any point of the manifold; a basis for an 
atlas is provided by the charts on 23 naturally induced from those on JVC by the representation of 
15 See, e.g., Penrose and Rindler (1984). 

16 In a more thorough treatment, one would characterize the way that the induction fixes a group action on the 
fibers, but we do not need to go into that for our purposes. Again, see Curiel (2009). 
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tensors on M as collections of components in M's coordinate systems. A cross-section is then the 
representation of a particular tensor field of that index structure on M. There is a natural induction 
in this case, i : Dm — > fixed by the pull-back action of a diffeomorphism <p of tensors on M. It is 
straightforward to show that i so defined is in fact an induction. Spinor bundles provide interesting 
examples of physically important bundles that cannot naturally be turned into geometric bundles. 

We are finally in a position to define concomitants. Let (551, M, 7Ti, t\) and (53 2 , M, 7r 2 , 12) be 
two bundles with the same base space. 17 

Definition 3.4 A mapping \ '■ ®i - * ^2 a concomitant if 



for all ui £ Hi and all £ Sjvt- 

In intuitive terms, a concomitant is a mapping between bundles that commutes with the action of the 
induced diffeomorphisms that lend the objects of the bundles their respective geometric characters. 
It is easy to see that \ must be fiber-preserving, in the sense that it maps fibers of 23 1 to fibers of 532. 
This captures the idea that the dependence of the one type of object on the other is strictly local, 
as the respecting of the actions of diffeomorphisms captures the idea that the mapping encodes an 
invariant relation. 

4 Jet Bundles and Higher-Order Concomitants 

Just as with ordinary functions from one Euclidean space to another, it seems plausible that the 
dependence encoded in a concomitant from one geometric bundle to another may take into account 
not only the value of the first geometrical structure at a point of the base space, but also "how 
that value is changing" in a neighborhood of that point, something like a generalized derivative of 
a geometrical structure on a manifold. The following construction is meant to capture in a precise 
sense the idea of such a generalized derivative in such a way so as to make it easy to generalize the 
idea of a concomitant to account for it. 

Fix a geometric fibre bundle (53, M, 7r, i), and the space of sections r[53]. Two sections 7,77 : 
M — > 53 osculate to first-order at p £ M if T7 and Tr\ agree in their action on T p M. 18 If (x l , v a ) are 
coordinates at the point -f(q) adapted to the bundle structure (as defined by the induction), then a 
coordinate representation of this relation is: 



for all i < dim(M) and a < dim(7r 1 [q\). This constitutes an equivalence relation on T[53]. A 1-jet 

17 One can generalize the definition of concomitants to cover the case of bundles over different base spaces, but we 
do not need this here. 

18 They osculate to zeroth-order at p if they map p to the same point in the domain. 
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with source q and target 7(g), written '^[7]', is such an equivalence class. The set of all 1-jets, 

= u 4 [7] 

gGM,7er[S] 

naturally inherits the structure of a diffcrcntiablc manifold. Let (0, U) be an adapted coordinate 
chart of 23 around 7(g), with the coordinate functions (V, /). Then the induced coordinate chart 
on J 1 23 is ((f) 1 , U 1 ) where 

^ = {j\[i\\l{q)€U} (4.1) 
and the coordinate functions associated with cf) 1 are (x l , v a , vf), where 

v a Vl[i]) = v a (7(q)) 



- - WW/ (-4 2) 

d(v a o 7) 



dx % 

where 7 is any member of ^[7]; this is well defined since all members of j 1 ^) agree on -f(q) and 
d(v a o 7) 



dx % 



by definition. 



g 

One can naturally fibre J lr B over M. The source projection a 1 : J lr B — > M, defined by 

gives J 1 23 the structure of a bundle space over the base space M, and in this case we write the 
bundle (J 1 23, M, cr 1 ). A section 7 of 23 naturally gives rise to a section ^[7] of J lr B, the first-order 
prolongation of that section: 

j 1 ^] : M -» [J #[ 7 ] 

such that tr 1 (j 1 [7]((7)) = g. (We assume for the sake of simplicity that global cross-sections exist; 
the emendations required to treat local cross-sections are trivial.) 

The points of J 1 23 may be thought of as coordinate- free representations of first-order Taylor 
expansions of sections of 23. To see this, consider the example of the trivial bundle (23, R 2 , tt) where 
23 = M 2 x M and n is projection onto the first factor. Fix global coordinates (a; 1 , x 2 , v 1 ) on 23, 
so that the induced (global) coordinates on J 1 23 are (a; 1 , a; 2 , v 1 , v 1 , v\). Then for any 1-jet jg[y], 
define the inhomogenous linear function 7 : M 2 — > R by 

lip) = ^ilip)) + v{(jt[l])iPl - 9l) + w 2(jg[7])(P2 - 92) 

where 7 G and p,g€l 2 with respective components (p\, P2) and (qi, q 2 ). Clearly 7 defines 

a cross-section of J 1 23 first-order osculant to 7 at p and so is a member of jg[y]; indeed, it is the 
unique globally defined, linear inhomogeneous map with this property. 
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A 2-jet is defined similarly, as an equivalence class of sections under the relation of having the 
same first and second partial-derivatives at a point. More precisely, 7,77 € r[25] osculate to second 
order at q € M if 7(g) = n{q) and 



9(v Q o 7) 



dx % 
d 2 {v a o 1 ) 



d(v a o rj) 



dx % dxi 



dx l 
d 2 {v a o V ) 



dx l dxi 



(4.3) 



One then defines J 2 25, et al, in the analogous ways. There is a natural projection from J 2 25 to J lr B, 
the truncation 6* 2,1 , characterized by "dropping the second-order terms in the Taylor expansion". In 
general, one has the natural truncation Q n > m : J n 25 — > J m 25 for all < m < n. 

An important fact for the present goal of defining concomitants is that the jet bundles of a 
geometric bundle are themselves naturally geometric bundles. Fix a geometric bundle (25, M, n, t) 
and a diffeomorphism (f> on M. Then l[<P] not only defines an action on points of 25, but, as a 
diffcomorphism itself on 25, it naturally defines an action on the cross-sections of 25 and thus on 
the 1-jets. The action can be characterized by the appropriate coordinate transformations of the 
formulae (4.2) and (4.3). It is easy to show that the mapping i 1 so specified from T)m to is 
an injective homomorphism and thus itself an induction; therefore, (J x 25, M, a , l 1 ) is a geometric 
fiber bundle. One defines inductions for higher-order jet bundles in the same way. 

We can now generalize our definition of concomitants. Let (25!, M, 7Ti, l) and (25 2 , M, 7r 2 , j) be 
two geometric fiber bundles over the manifold M. 

Definition 4.1 A zeroth-order concomitant from 25i to 25 2 is a smooth mapping \ '■ 23 1 — > 25 2 such 
that 

(Vp e Si)(V^ e 21 m ){j(0)(x(p)) = xW0)(p))} 

Definition 4.2 y4n n th -order concomitant > 0) from 25i to 25 2 is a smooth mapping \ '■ J n< S>i — * 
25 2 such that 

1. (Vu g J"25!)(V0 e a M ) = X(* n (0)(«)) 

^. i/iere is no (n — l) th concomitant \' : J n ~ lr i>\ — » 25 2 satisfying 

(VueJ n V 1 )x(u) = x'(O n ' n - 1 (u)) 

I will sometimes refer to a zeroth-order concomitant simply as a concomitant for the sake of brevity. 
By another convenient abuse of terminology, I will often refer to the range of the concomitant 
mapping itself as 'the concomitant' of the domain. It will be of physical interest in §6 to consider 
the way that concomitants interact with multiplication by a scalar field. (Since we consider in this 
paper only concomitants of linear and affine objects, multiplication of the object by a scalar field is 
always defined.) In particular, let us say that a concomitant is homogeneous of weight w if for any 
constant scalar field a 

An important property of concomitants is that, in a limited sense, they are transitive. 
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Proposition 4.3 If Xi '■ J nr &i — * 2>2 is an n th -order concomitant and \2 ■ — » ^3 *s a smooth 
mapping, where 23 1; H2 an d 23 3 are geometric bundles over the same base space, then \2 Xi * s arl 
n th -order concomitant if and only if X2 is a zeroth-order concomitant. 

This follows easily from the fact that inductions are injective homomorphisms and concomitants 
respect the fibers. 

5 Concomitants of the Metric 

As a specific example that will be of use in what follows, consider the geometric fiber bundle 
(23g, M, 7Tg, Lg), with M a 4-dimensional, paracompact, connected, smooth manifold (i.e., a can- 
didate spacetimc manifold), the fibers of CBg diffcomorphic to the set of Lorcntz metrics at each 
point of JVC, all of the same signature (+, — , — , — ), and tg the induction defined by the natural 
pull-back. Since the set of Lorentz metrics in the tangent plane over a point of a 4-dimensional 
manifold, all of the same signature, is a 10-dimensional topological space, 19 the bundle space 23g is 
a 14-dimcnsional manifold. A cross-section of this bundle is equivalent to a Lorcntz metric field on 
the manifold. 

The following proposition precisely captures the statement one sometimes hears that there is no 
scalar or tensorial quantity one can form depending only on the metric and its first-order partial 
derivatives at a point of a manifold. 

Proposition 5.1 There is no first-order concomitant from "Bg to any tensor bundle over M. 

To prove this, it suffices to remark that, given any spacetime (M, g a b) and any two points p,p' € M, 
there are coordinate neighborhoods U of p and U 1 of p 1 and a diffcomorphism <f> : JA, — ► M, such that 
0(p) = p' 1 4^{9'ab) = 9°.b at P, and 4>$(d' a gbc) = d a gbc at p, where d a (d' a ) is the ordinary derivative 
operator associated with the coordinate system on U (If), and <ffi is the map naturally induced by 
the pull-back action of (f>. 

This is not to say, however, that no information of interest is contained in J lr E>g- Indeed, two 
metrics g a b and h a b are first-order osculant at a point if and only if they have the same associated 
derivative operator at that point. To see this, first note that, if they osculate to first order at that 
point, then V a (9bc ~ ^&c) = at that point for all derivative operators. Thus, for the derivative 
operator V a associated with, say, g a b, ^a(gtc — hb c ) = 0, but V a gbc = 0, so V a hb c = at that point 
as well. Similarly, if the two metrics are equal and share the same associated derivative operator V a 
at a point, then V a (gb c — hb c ) = at that point for all derivative operators, since their difference 
will be identically annihilated by V , and g a b = h a b at the point by assumption. Thus they are first- 
order osculant at that point and so in the same 1-jet. This proves that all and only geometrically 
relevant information contained in the 1-jets of Lorentz metrics on M is encoded in the fiber bundle 
over spacetime the values of the fibers of which are ordered pairs consisting of a metric and the 
metric's associated derivative operator at a spacetime point. 

19 In fact, it is homeomorphic to a connected, convex, open subset — an open cone with vertex at the origin — in R 10 . 
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The second jet bundle over Sg has a similarly interesting structure. Clearly, if two metrics are 
in the same 2-jet, then they have the same Riemann tensor at the point associated with the 2-jet, 
since they have the same partial-derivatives up to second order at the point. Assume now that two 
metrics are in the same 1-jet and have the same Riemann tensor at the associated spacetime point. 
If it follows that they are in the same 2-jet, then essentially all and only geometrically relevant 
information contained in the 2-jets of Lorentz metrics on M is encoded in the fiber bundle over 
spacetime the points of the fibers of which are ordered triplets consisting of a metric, the metric's 
associated derivative operator and the metric's Riemann tensor at a spacetime point. To demonstrate 
this, it suffices to show that if two Levi-Civita connections agree on their respective Riemann tensors 
at a point, then the two associated derivative operators are in the same 1-jet of the bundle whose 
base-space is M and whose fibers consist of the affine spaces of derivative operators at the points of 
M (because they will then agree on the values of first-partial derivatives of their Christoffel symbols 
at that point in any coordinate system as well as agreeing in the values of the Christoffel symbols 
themselves, and thus will be in the 2-jet of the same metric at that point). 

Assume that, at a point p of spacetime, g ab — g ab , V a = V a (the respective derivative operators), 
and R a bcd = R a bcd (the respective Riemann tensors). Let C a bc be the symmetric diffcrcncc-tensor 
between V a and V , which is itself at p by assumption. Then by definition V[ b V c ]£ a = R a bcnS, n 
for any vector £°, and so at p by assumption 

= v a (v fc r + c c fc „C)-v b v a £ c 

= V a V 6 £ c + v a (c c fc „o - v 6 v„e c - c c fc „v a c + C" ba V„£ c 

but V a V b £ c — V b V a £ c = R a bcn£, n and C a bc = 0, so expanding the only remaining term gives 

rv a c c b „ = o 

for arbitrary £ a and thus V a C fc c d = at p; by the analogous computation, V a C b c d = as well. It 
follows immediately that V a and V a are in the same 1-jet of affine bundles over p. We have proved 

Theorem 5.2 J lr Bg is naturally diffeomorphic to the geometric fiber bundle over M whose fibers 
consist of pairs {g a b, V a ), where g a b is the value of a Lorentz metric field at a point ofM, and V a is 
the value of the covariant derivative operator associated with g a b at that point, the induction being 
defined by the natural pull-back. J 2l Bg is naturally diffeomorphic to the geometric fiber bundle over 
M whose fibers consist of triplets (g a b, V a , R a bc d ), where g a b is the value of a Lorentz metric field at 
a point ofJA, and V a and R a bc d are respectively the covariant derivative operator and the Riemann 
tensor associated with g a b at that point, the induction being defined by the natural pull-back. 

It follows immediately that there is a first-order concomitant from CBg to the geometric bundle 
(23y, M, 7Tv, iv) °f derivative operators, viz., the mapping that takes each Lorentz metric to its 
associated derivative operator; likewise, there is a second-order concomitant from 'Bg to the geometric 
bundle (35Rj em , M, ^ja em , £Riem) of tensors with the same index structure and symmetries as the 
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Ricmann tensor, viz., the mapping that takes each Lorentz metric to its associated Riemann tensor. 
This is the precise sense in which the Riemann tensor associated with a given Lorentz metric is "a 
function of the metric and its partial derivatives up to second order" . 

6 No Gravitational Stress-Energy Tensor Exists 

We are finally in a position to prove the main result of the paper, the nonexistence of a gravitational 
stress-energy tensor. In order to formulate and prove a result having that proposition as its natural 
interpretation, one must first lay down some natural conditions on the proposed object, to show that 
no such object exists satisfying the conditions. In general relativity, the invariant representation of 
energetic quantities is always in the form of a stress-energy tensor, viz., a two-index, symmetric, 
divergence-free tensor. Not just any such tensor will do, however, for that gives only the baldest of 
formal characterizations of it. From a physical point of view, at a minimum the object must have 
the dimension of stress-energy for it to count as a stress-energy tensor. That it have the dimension 
of stress-energy is what allows one to add two of them together in a physically meaningful way to 
derive the physical sum of total stress-energy from two different sources. In classical mechanics, for 
instance, both velocity and spatial position have the form of a three-dimensional vector, and so their 
formal sum is well defined, but it makes no physical sense to add a velocity to a position because 
the one has dimension of length/time and the other the dimension of length. 

An essential, defining characteristic of energy in classical physics is its obeying some formulation 
of the First Law of Thermodynamics. The formulation of the First Law I rely on is somewhat 
unorthodox: that all forms of stress-energy are in principle ultimately fungible — any form of energy 
can in principle be transformed into any other form — not necessarily that there is some absolute 
measure of the total energy contained in a system or set of systems that is constant over time. 
In more precise terms, this means that all forms of energy must be represented by mathematical 
structures that allow one to define appropriate operations of addition and subtraction among them, 
which the canonical form of the stress-energy does allow for. 20 I prefer this formulation of the First 
Law in general relativity because there will not be in a general cosmological context any well-defined 
global energetic quantity that one can try to formulate a conservation principle for. 21 In so far as 
one wants to hold on to some principle like the classical First Law in a relativistic context, therefore, 
I see no other way of doing it besides formulating it in terms of fungibility. (If one likes, one can 
take the fungibility condition as a necessary criterion for any more traditional conservation law.) 
This idea is what the demand that all stress-energy tensors, no matter the source, have the same 
physical dimension is intended to capture. 22 

20 Note that this is a requirement even if one takes a more traditional view of the First Law as making a statement 
about global conservation of an absolute magnitude measuring a physical quantity. 

21 It is interesting to note that, for similar reasons, there will in general be no available global statement of the 
relativistic analogue of the Second Law either. 

22 For what it's worth, this conception has strong historical warrant — Einstein (implicitly) used a very similar 
formulation in one of his first papers laying out and justifying the general theory (Einstein 1916, p. 149): 

It must be admitted that this introduction of the energy-tensor of matter is not justified by the relativity 
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To sum up, the stress-energy tensor encodes in general relativity all there is to know of ponderable 
energetic phenomena at a spacetime point: 

1. it has 10 independent components representing with respect to a fixed pseudo-orthonormal 
frame, say, the 6 components of the classical stress-tensor, the 3 components of linear momen- 
tum and the scalar energy density of the ponderable field at that point 

2. that it has two covariant indices represents the fact that it defines a linear mapping from 
timelike vectors at the point ("worldline of an observer") to covectors at that point ("4- 
momentum covector of the field as measured by that observer"), and so defines a bi- linear 
mapping from pairs of timelike vectors to a scalar density at that point ( "scalar energy density 
of the field as measured by that observer"), because energetic phenomena, crudely speaking, 
are marked by the fact that they are quadratic in velocity and momental phenomena by being 
linear in velocity 

3. that it is symmetric represents, "in the limit of the infinitesimal", the classical principle of 
the conservation of angular momentum and the special relativistic equivalence of momentum- 
density flux and scalar energy density 

4. that it is covariantly divergence-free represents the fact that, "in the limit of the infinitesimal" , 
the classical principles of energy and linear momentum conservation are obeyed 

5. the localization of ponderable stress-energy is embodied in the fact that the object representing 
it is a tensor, a multi-linear map acting only on the tangent plane of the point it is associated 
with 

6. finally, the thermodynamic fungibility of energetic phenomena is represented by the fact that 
the set of stress-energy tensors forms a vector space — the sum and difference of any two is 
itself a possible stress-energy tensor — all having the same physical dimension 

Consequently, the appropriate mathematical representation of localized gravitational stress-energy, 
if there is one, is a two covariant-index, symmetric, covariantly divergence-free tensor having the 
physical dimension of stress-energy. 

In order to make precise the idea of having the physical dimension of stress-energy, recall that 
in general relativity all the fundamental units one uses to define stress-energy, namely time, length 
and mass, can themselves be defined using only the unit of time; these are so-called geometrized 
units. For time, this is trivially true: stipulate, say, that a time-unit is the time it takes a certain 
kind of atom to vibrate a certain number of times under certain conditions. A unit of length is then 
defined as that in which light travels in vacuo in one time-unit. A unit of mass is defined as that of 

postulate alone. For this reason we have here deduced it from the requirement that the energy of the 
gravitational field shall act gravitatively in the same way as any other kind of energy. 

Moller (1962) also stresses the fact that the formulation of integral conservation laws in general relativity based on 
pscudo-tensorial quantities depends crucially on the assumption that gravitational energy, such as it is, shares as 
many properties as possible with the energy of ponderable matter. 
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which two, placed one length- unit apart, will induce in each other by dint of their mutual gravitation 
alone an acceleration towards each other of one length-unit per time-unit per time-unit. 23 These 
definitions of the units of mass and length guarantee that they scale in precisely the same manner 
as the time-unit when new units of time are chosen by multiplying the time-unit by some fixed 
real number A - 5. (The reason for the inverse square-root will become clear in a moment). Thus, 
a duration of t time-units would become t\~^ of the new units; an interval of d units of length 
would likewise become <iA~2 in the new units, and m units of mass would become mA"? of the new 
units. This justifies treating all three of these units as "the same" , and so expressing acceleration, 
say, in inverse time-units. To multiply the length of all timelike vectors representing an interval 
of time by A~5 ; however, is equivalent to multiplying the metric by A (and so the inverse metric 
by A -1 ), and indeed such a multiplication is the standard way one represents a change of units in 
general relativity. This makes physical sense as the way to capture the idea of physical dimension: 
all physical units, the ones composing the dimension of any physical quantity, are geometrized in 
general relativity, in the most natural formulation, and so depend only on the scale of the metric 
itself. 

Now, the proper dimension of a stress-energy tensor can be determined by the demand that the 

Einstein field equation, G a b = jT a b, where 7 is Newton's gravitational constant, remain satisfied 

when one rescales the metric by a constant factor. 7 has dimension (length) ^ gQ 

(mass) (time) 

geometrized units does not change under a constant rescaling of the metric. Thus T a b ought to 
transform exactly as G a b under a constant rescaling of the metric. A simple calculation shows that 
G a b = Rab — \Rgab remains unchanged under such a rescaling. Thus, a necessary condition for 
a tensor to represent stress-energy is that it remain unchanged under a constant rescaling of the 
metric. It follows that the concomitant at issue must be homogeneous of degree in the metric. 

It is clear that a gravitational stress-energy tensor, if there is such a thing, ought to be a 
concomitant of the metric, but the order of the required concomitant is not a priori obvious. In 
fact, the way a homogeneous concomitant of the metric transforms when the metric is multiplied 
by a constant factor suffices to fix the differential order of that concomitant. 24 This can be seen as 
follows, as exemplified by the case of a two covariant-index, homogeneous concomitant S a b of the 
metric. A simple calculation based on definition 4.2 and on the fact that the concomitant must be 
homogeneous shows that the value at a point p e M of an n th -order concomitant S a b can be written 

23 This definition may appear circular, in that it would seem to require a unit of mass in the first place before one 
could say that bodies were of the same mass. I think the circularity can be mitigated by using two bodies for which 
there are strong prior grounds for positing that they are of equal mass, e.g., two fundamental particles of the same 
type. It also suffers from a fundamental lack of rigor that the definition of length does not suffer from. In order to 
make the definition rigorous, one would have to show, e.g., that there exists a solution of the Einstein Field-Equation 
(approximately) representing two particles in otherwise empty space (as defined by the form of T ab ) — viz., two timelike 
geodesies — such that, if on a spacelike hypersurface at which they both intersect 1 unit of length apart (as defined 
on the hypersurface with respect to cither) they accelerate towards each other (as defined by relative acceleration of 
the geodesies) one unit length per unit time squared, then the product of the masses of the particles is 1. I will just 
assume, for the purposes of this paper, that such solutions exist. 

24 1 thank Robert Geroch for pointing this out to me. 
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in the general form 

S ab = £ k a g«* ...g xr (Vi"i V) . . . (v< n 'W) (6.1) 

a 

where: V a is any derivative operator at p other than the one naturally associated with g a b\ the x's 
arc dummy abstract indices; < vi™ , ' ) ' stands for rij iterations of that derivative operator (obviously 
each with a different abstract index); a takes its values in the set of all permutations of all sets of 
positive integers {m, . . . , n,} that sum to n, so i can range in value from 1 to n; the exponents of 
the derivative operators in each summand themselves take their values from a, i.e., they are such 
that n\ + • • • + rii = n; there is exactly one summand for which ri\ — n (which makes it an n th -order 
concomitant); for each a, k a is a constant; and there are just enough of the inverse metrics in each 
summand to contract all the covariant indices but a and b. 
Now, a combinatorial calculation shows the following: 

Proposition 6.1 If, for n>2, S a b is an n th -order homogeneous concomitant of g a b, then to rescale 
the metric by the constant real number A multiplies S a b by A"~ 2 . 

In other words, the only such homogeneous n th -order concomitants must be of weight A — 2. So if one 
knew that S a b were multiplied by, say, A 4 when the metric was rescaled by A, one would know that it 
had to be a sixth-order concomitant. In particular, S a b does not rescale when g a b — > \g a b if and only 
if it is a second-order concomitant of g a b- 25 Thus, such a tensor has the physical dimension of stress- 
energy if and only if it is a second-order concomitant. In consequence, the only viable candidates for 
gravitational stress-energy tensors are two covariant-index, symmetric, divergence- free, second-order 
concomitants of the metric. 

We can now state and prove the main result. 

Theorem 6.2 The only two covariant-index, symmetric, divergence-free, second-order concomitants 
of the metric are constant multiples of the Einstein tensor. 

(Note that homogeneity of the concomitant does not need to be assumed; it was used only to fix the 
required order of the concomitant at issue.) Before proving the theorem, I remark that it does bear 
the required natural interpretation, for the Einstein tensor is not an appropriate candidate for the 
representation of gravitational stress-energy: the Einstein tensor will be zero in a spacetime having 
vanishing Ricci tensor but non-trivial Weyl tensor (i.e., one that has only non-trivial conformal 
curvature); such spacetimes, however, can manifest phenomena, e.g., gravitational radiation, that 

25 Note that the exponent (n — 2) in this result depends crucially on the fact that S a i, has only two indices, both 
covariant. This result generalizes as follows: if S a ' i is an n th -order concomitant of g ab such that the number of 
its covariant indices minus the number of its contravariant indices is 2, then, under constant rescaling of the metric 
by the real number A, S a will get multiplied by a factor of A n ~ 2 . Similar propositions, with different exponents 
resulting for A, can be proven in the same way for tensor concomitants of the metric of any index structure. A slight 
variation of the argument, moreover, shows that there does not in exist in general a homogeneous concomitant of a 
given order from a tensor of a given index structure to one of another structure — one may not be able to get the 
number and type of the indices right by contraction and tensor multiplication alone. 
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one naturally wants to say possess gravitational energy in some (necessarily non-localized) form or 
other. 26 

Now, to prove the theorem, note first that it follows from theorem 5.2 and proposition 4.3 that 
any candidate gravitational stress-energy tensor must be a zeroth-order concomitant of "B]^j em , the 
geometric bundle of Riemann tensors over spacetime. (One can take this as a precise statement of 
the fact that any gravitational stress-energy tensor ought to "depend on the curvature" , as I argued 
in §2.) 

One can form from algebraic combinations of the Riemann tensor and the metric 14 independent 
scalar invariants, 27 which manifestly are themselves concomitants of the metric of order at least 
second. The important point about these scalar invariants for this argument is that the result 
of taking the variation of any of them with respect to the metric yields a two covariant-index, 
symmetric, divergence-free tensor that is itself a concomitant of the metric of the same order as the 
scalar invariant; moreover, all two covariant-index, symmetric, divergence-free concomitants of the 
metric, of all orders second and higher, can be derived in this fashion. 

To make this claim precise, let 6 g be the space of scalar invariants of R a bc d on M with associated 
metric g a b, and g(X) a b be a smooth one-parameter family of metrics defined for A e (—1, 1) such that 
g{0)ab = 9ab- I will be interested in the inverse metric g ab and the associated smooth one-parameter 
family g{X) ab ; denote by l Sg ab \ Define the functional S : 6 g -» M by 

S[T,} = / Sv 7 ^? e abcd 
Jm 

for E e 6g, where e abc d is a fixed 4-form (not necessarily the one associated with g a b) and \J~^g 
is the determinant of g ab in a coordinate system in which the only components of e a bcd are 0, 1 or 
— I. 28 Clearly | A _ exists for all smooth one-parameter families g(X) a b', denote it by l 5S\ Then 
the variation of E with respect to g ab is the unique symmetric tensor field such that 

5S= S nm 6g nm y/=g e abcd 

for all smooth one-parameter families g(\) a b- A straightforward calculation shows that the variation 
of any scalar curvature invariant has vanishing covariant divergence with respect to the derivative op- 
erator associated with the Riemann tensor given metric; furthermore, the variation is a concomitant 
of the associated metric of the same order as the scalar curvature invariant itself. 

The tensor fields resulting from varying all such scalar curvature invariants with respect to the 
fixed metric form a real vector space over each point of spacetime of dimension 6, the same as that 
of all two covariant-index, symmetric, tensors divergence-free with respect to the fixed metric. From 



26 As an historical aside, it is interesting to note that early in the debate on gravitational energy in general relativity 
Lorentz (1916) and Levi-Civita (1917) proposed that the Einstein tensor be thought of as the gravitational stress- 
energy tensor. Einstein criticized the proposal on the grounds that this would result in attributing zero total energy 
to any closed system. 

27 Cf. Haskins (1902). 

28 If M is not compact, or if the integral does not converge for some other reason, then it ought to be taken over 
some appropriate compact subset of JYt, on the boundaries of which 5g ab is to be set to zero. 
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this, in conjunction with the fact that S is a linear functional, it follows that all two covariant-indcx, 
symmetric, divergence- free concomitants of any order of the Riemann tensor (and thus of the metric 
of any order second or higher) can be formed by taking the variation of all scalar curvature invariants 
of the Riemann tensor with respect to the metric. Thus any candidate gravitational stress-energy 
tensor must be the variation with respect to the metric of a scalar curvature invariant. 

According to the definition of a homogeneous concomitant, a scalar curvature invariant that is, 
e.g., quadratic in the Riemann tensor, that is to say, quadratic in second derivatives of the metric, 
is a fourth-order concomitant of the metric. The rule is that the order of the concomitant is the 
sum of the exponents of the derivative operators when the concomitant is represented in the form 
of equation (6.1). 29 Consequently a scalar curvature invariant is a second-order concomitant of the 
metric — and so its variation with respect to the metric yields a second-order concomitant of the 
metric — if and only if it is linear in the Riemann tensor. The only scalar curvature invariants that 
are linear in the Riemann tensor are constant multiples of the Gaussian scalar curvature itself, R, 
viz., the result of contracting the two indices of the Ricci tensor. The tensor that results from taking 
the variation of R with respect to the metric is the Einstein tensor, which proves the theorem. 
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